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A  comprehensive  mathematical  model  for  the  planar  solid  oxide  fuel  cell  (SOFC)  stack  is  presented.  The 
model  couples  the  intricate  interdependency  among  the  ionic  conduction,  electronic  conduction,  gas 
transport  and  the  electrochemical  reaction  and  takes  into  account  the  contact  resistance  between  the 
electrode  and  interconnect  rib  and  the  dependence  of  the  effective  electrode  properties  on  the  microstruc¬ 
ture  parameters  of  the  porous  electrode.  The  validity  of  the  mathematical  model  is  demonstrated  by  the 
excellent  agreement  between  the  numerical  and  experimental  I-V  curves.  Based  on  a  standard  set  of 
model  parameters,  the  cell  performance  is  examined  by  systematically  varying  the  rib  width,  contact 
resistance,  fuel  composition  electrode  properties  and  the  pitch  width.  The  results  show  conclusively  that 
the  cell  output  depends  strongly  on  the  rib  width  and  a  suitable  choice  of  the  rib  width  is  important  for 
realizing  the  potential  of  a  SOFC  stack.  The  optimal  rib  width  is  only  sensitive  to  the  contact  resistance 
and  the  pitch  width,  but  the  optimal  results  for  the  anode  and  cathode  ribs  are  quite  different.  Both  the 
optimal  anode  and  cathode  rib  widths  depend  linearly  on  the  contact  resistance  and  the  pitch  width  and 
the  parameters  for  the  linearity  are  given. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  anode-supported  planar  solid  oxide  fuel  cell  (pSOFC)  is  a 
promising  new  energy  technology  and  has  attracted  increasing 
attention  in  the  past  years.  It  has  long  been  reported  that  a  single 
cell  of  pSOFC  may  produce  a  power  density  of  1  Wcnrr2  or  higher 
at  the  operating  temperature  of  700-800  °C  [1-6].  Although  the 
performance  of  pSOFC  at  the  single  cell  level  is  exciting,  multiple 
cells  must  be  connected  in  series  by  interconnect  material  to  form 
a  fuel  cell  stack  for  practical  applications.  Unfortunately,  the  per¬ 
formance  of  a  pSOFC  cell  at  a  working  stack  level  rarely  exceeds 
a  third  of  1  Wcnrr2  [7].  The  cell  performance  for  small  laboratory 
stacks  can  be  higher,  but  significant  reductions  of  power  density 
due  to  the  stack  assembly  are  widely  observed  [8]. 

There  are  various  factors  that  may  affect  the  stack  cell  perfor¬ 
mance.  One  of  the  primary  causes  accounting  for  the  lowered  cell 
performance  at  the  stack  level  is  the  additional  losses  caused  by 
the  interconnect  geometry  [8,9].  Small  grooves  in  interconnect  are 
commonly  used  as  gas  channels  to  carry  the  fuel  and  air  flows.  The 
ribs,  which  separate  and  define  the  channels,  make  direct  contact 
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with  the  electrodes.  On  the  technical  side,  the  channels  and  the 
ribs  cannot  be  too  small  and  the  contact  resistance  between  the 
ribs  and  the  electrodes  is  part  of  the  intrinsic  property  of  the  inter¬ 
connect  material  [10].  In  designing  the  layer  architecture,  there  is 
a  tradeoff  that  must  be  considered  between  the  rib  and  channel 
sizes.  On  one  hand,  ribs  covering  bigger  fraction  of  the  cell  area 
may  reduce  the  interface  resistance  to  current  flow  by  increasing 
the  electrode-interconnect  contact  area.  Such  ribs  will  give  a  bet¬ 
ter  conduction  of  the  electrical  current  and  reduce  ohmic  losses. 
On  the  other  hand,  the  chemical  species  do  not  diffuse  as  well 
underneath  wide  ribs.  Narrow  ribs  and  ribs  covering  small  fraction 
of  the  cell  area  are  needed  to  facilitate  more  uniform  distribu¬ 
tion  of  reactive  gases  across  the  area  of  the  electrolyte  surface 
and  thus  to  promote  electrochemical  performance.  The  implica¬ 
tion  of  the  tradeoff  to  the  cell  performance  can  be  very  significant. 
It  was  reported  that  the  power  density  was  0.76  W  cm-2  at  0.7  V 
operating  voltage  for  the  combination  of  2.6  mm  channels  and 
2.4  mm  ribs,  when  that  for  4  mm  channels  and  1  mm  ribs  was  as 
high  as  1.03  W cm-2  [8].  Clearly,  how  to  configure  a  channel  net¬ 
work  to  improve  the  stack  cell  performance  is  highly  consequential 
for  the  practical  application  [11]  and  is  an  important  task  of  the 
pSOFC  designer. 

There  have  been  a  few  theoretical  studies  investigating  the 
impact  of  the  rib  size  on  the  performance  of  solid  oxide  fuel  cells. 
Ferguson  et  al.  [12]  gave  a  numerical  example  for  the  rib  size  effect 
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on  the  performance  of  a  specific  electrolyte-supported  SOFC.  Tan¬ 
ner  and  Virkar  proposed  a  simplified  model  to  examine  the  effect 
of  the  interconnect  design  on  the  effective  Ohmic  resistance  of 
the  membrane-electrode  assembly  in  an  SOFC  stack  [13].  Lin  et  al. 

[14]  described  a  phenomenological  model  and  provided  analyt¬ 
ical  expressions  to  estimate  the  rib  effects  on  the  concentration 
and  ohimc  polarizations  of  anode-supported  SOFC  stacks.  Ji  et  al. 

[15]  showed  that  the  terminal  output  of  a  SOFC  stack  depended 
strongly  on  the  contact  resistance.  Jeon  et  al.  [16]  described  a 
microstructural  model  and  examined  the  influence  of  the  rib  and 
pitch  widths  and  the  electrode-interconnect  contact  area  specific 
resistance  ( ASRcontact )  on  the  stack-cell  performance.  Based  on  sys¬ 
tematic  examination  of  the  stack  cell  performance  as  a  function  of 
various  influencing  factors,  Liu  et  al.  [  1 7 ]  proposed  a  simple  expres¬ 
sion  for  the  optimal  rib  width  as  a  function  of  the  pitch  width  and 
ASRcontact  that  was  independent  of  the  porosity,  layer  thickness  and 
conductivity  of  the  electrode  and  easy  to  use  for  the  SOFC  engineers. 

Unfortunately,  the  existing  numerical  studies  on  the  effect  of  the 
rib  size  on  the  stack  cell  performance  have  been  carried  out  with  the 
assumption  of  equal  width  of  the  anode  and  cathode  ribs  and  equal 
ASRcontact  for  the  anode-interconnect  and  cathode-interconnect 
interfaces.  As  a  result,  it  is  difficult  to  distinguish  between  the 
impacts  of  the  anode  rib  and  the  cathode  rib.  As  demonstrated 
experimentally  by  Kornely  et  al.  [8],  however,  ASRcontact  for  the 
anode-interconnect  and  cathode-interconnect  interfaces  can  be 
very  different  and  the  impacts  of  the  anode  rib  and  the  cathode  rib 
on  the  stack  cell  performance  are  also  very  different.  In  fact,  previ¬ 
ous  numerical  studies  have  already  shown  very  different  effects  of 
the  cathode  and  anode  ribs  on  the  mass  transport  in  the  electrodes. 
For  example,  Liu  et  al.  [17]  reported  that  the  anode  gas  transport 
of  an  anode-supported  SOFC  was  only  mildly  affected  by  the  anode 
rib,  while  an  oxygen  depletion  zone  of  0.46  mm  was  found  with  a 
cathode  rib  width  of  only  0.8  mm.  The  existing  experimental  and 
theoretical  findings  indicate  clearly  that  the  optimal  design  of  the 
stack  cell  should  use  different  widths  for  the  anode  and  the  cathode 
ribs. 

Building  upon  the  foundation  of  significant  earlier  literatures 
[18,19],  this  paper  describes  a  two-dimensional  microstructure 
based  multi-physics  model  for  SOFCs  with  composite  electrodes. 
The  model  couples  the  intricate  interdependency  among  the  ionic 
conduction,  electronic  conduction,  gas  transport  in  porous  elec¬ 
trodes  and  the  electrochemical  reaction  processes  in  the  functional 
layers  and  at  the  electrode/electrolyte  interfaces.  The  impacts  of 
the  anode  and  cathode  ribs  are  examined  independently  and  for 
different  electrode-interconnect  contact  resistances,  fuel  compo¬ 
sitions,  electrode  porosities,  electrode  thicknesses  and  electrode 
conductivities.  The  optimal  cathode  and  anode  rib  sizes  for  different 
ASRcontact  are  obtained  and  parameterized  to  provide  an  easy-to-use 
guidance  for  designing  the  rib-channel  layout. 


2.  Method 

2.2.  Physical  model 

A  repeating  cell  unit  of  a  pSOFC  stack  is  shown  schemati¬ 
cally  in  Fig.  la.  The  repeating  unit  consists  of  (i)  cathode-side 
interconnect  plate  and  the  air  channels,  (ii)  positive  electrode 
(cathode)-electrolyte-negative  electrode  (anode)  assembly,  or 
PEN  for  short,  (iii)  anode-side  interconnect  plate  and  the  fuel  chan¬ 
nels.  A  PEN  assembly  with  composite  electrodes  often  consists 
of  five  layers:  (1)  an  anode  support  layer  (ASL)  that  also  acts  as 
the  anode  current  collector,  (2)  an  anode  functional  layer  (AFL) 
for  reducing  the  activation  losses  by  affording  abundant  three 
phase  boundaries  (TPBs),  (3)  an  electrolyte  layer,  (4)  a  cathode 
functional  layer  (CFL),  and  (5)  a  cathode  current  collector  layer 


Fig.  1.  Schematic  of  (a)  an  anode-supported  planar  SOFC  stack  cell  unit,  (b)  a  cross- 
section  of  the  SOFC  stack  cell  unit. 

(CCCL).  A  cross  section  of  the  repeating  stack  cell  unit  is  shown  in 
Fig.  lb. 

For  simplicity,  the  stack  cell  is  assumed  to  be  in  an  isothermal 
steady  state.  As  demonstrated  before  [17,19],  a  2D  model  is  equiv¬ 
alent  to  a  3D  model  when  discussing  the  effects  of  rib  sizes  on  the 
fuel  cell  performance  as  they  provide  essentially  the  same  results. 
Therefore,  a  2D  model  (Fig.  lb)  is  used  here.  Moreover,  we  assume 
the  electrical  conductivity  of  the  interconnect  is  sufficiently  high 
and  the  ohmic  loss  inside  the  interconnect  (excluding  the  ohmic 
loss  caused  by  the  electrode-interconnect  contact  resistance)  is 
negligible.  Taking  into  account  the  equivalency  of  the  repeating 
rib-channel  pairs  (pitches),  the  computational  domain  consists  of 
only  half  of  the  repeating  rib-channel  pair,  shown  in  Fig.  lb  as  inside 
the  red  square.  The  square  computational  domain  has  a  width  of 
dpitch  (the  sum  of  one  half  of  the  rib  width,  d rib,  and  one  half  of 
the  channel  width,  dchanne/ )  and  a  height  of  the  total  PEN  thick¬ 
ness.  The  computational  domain  represents  a  repeating  unit  in  a 
pSOFC  cell  of  any  size.  Notice  that  we  have  implicitly  assumed  here 
that  the  anode  and  cathode  pitch  widths  are  equal,  even  though 
the  widths  for  the  anode  rib  (cPJ )  and  cathode  rib  (cFab)  may  be 
different.  As  the  gas  transports  inside  the  anode  and  the  cathode 
are  practically  independent  from  each  other,  the  above  implicit 
assumption  for  numerical  convenience  should  not  put  any  practi¬ 
cal  limitation  on  the  modeling  result.  Indeed,  previous  3D  modeling 
studies  [19]  have  shown  that  the  rib  effects  on  the  overall  stack  cell 
output  are  basically  the  same  whether  the  stack  cell  has  a  co-flow, 
counter-flow  or  cross-flow  configuration. 

The  actual  cell  operating  voltage  ( Vop )  is  lower  than  the  ideal 
cell  potential  (Nernst  potential)  due  to  various  polarizations  such 
as  activation  polarization,  ohmic  polarization  and  concentration 
polarization.  For  the  physical  model  described  above,  the  Vop  of 
the  stack  cell  may  be  formally  expressed  as: 

\r  c  y,an  ^ an  7nan  ^ an  ^el  7nca 

op  —  0  'Iasr  'leone  'I  act  'l0hm  'l0hm  'l0hm 

-Ct-^nc-^R  (1) 

where  F0  is  the  Nernst  potential,  ij™R  (^R)  the  anode-rib  (cathode- 
rib)  interface  overpotential  due  to  the  contact  resistance  at  the 
material  boundary,  r\™nc  (ijcc“nc)  the  anode  (cathode)  concentra¬ 
tion  polarization  due  to  the  resistance  to  the  transport  of  fuel 
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(air)  species  through  the  porous  electrode,  (r)caact)  the  anode 
(cathode)  activation  polarization  due  to  the  energy  barrier  for 
the  electrochemical  reactions,  rf^  ,  rjeolhm  and  qcahm  the  ohmic 
polarizations  in  the  anode,  electrolyte  and  cathode  associated 
with  the  ohmic  resistances  of  the  electrolyte  and  electrodes.  The 
detailed  model  for  each  of  the  polarizations  will  be  described 
later. 


2.2.  Governing  equations  for  gas  transport  in  porous  electrodes 


The  mass  transport  of  gas  species  in  a  porous  medium  is  complex 
and  includes  three  distinct  mass-transfer  mechanisms:  Knudsen 
diffusion,  molecular  diffusion,  as  well  as  viscous  flow.  Thus,  the 
dusty  gas  model  is  required  to  estimate  accurately  the  gas  transport 
in  SOFC  porous  electrodes  [20-22].  For  a  binary  component  system 
(fuel  with  H2  and  H20,  or  air  with  02  and  N2),  the  total  molar  flux 
of  species  i  (i  =  1,  2)  may  be  written  as  [16,22,23], 


Nj  = 


neff  neff 
D\2  Dil( 


Dil  +xlDfK+x2D 


12 


eff  VC|' 
IK 


U\KU2K 


RTctot(De/f  +  xiDg  +  x2DfK)  M 


Vp 


2/C 


1/C 


—D  Vc-  —  C  ^  —  Ndiffusion  _|_  ^convection 

fi  1 


(2) 


where  c,-  is  the  molar  concentration  of  species  i,  ctot  (=ci+c2) 
the  total  molar  concentration  of  the  binary  mixture,  xI(=cI7ctot) 
the  molar  fraction  of  species  /,  R  the  universal  gas  constant,  T 
the  absolute  temperature,  p  the  total  gas  pressure,  k  the  perme¬ 
ability  coefficient,  /x  the  viscosity  coefficient,  the  effective 

Knudsen  diffusion  coefficient  of  species  i  and  the  effective 

binary  diffusion  coefficient.  D,(=  (D^D^)/(D^  +x\ +  x2D^)) 
is  the  equivalent  diffusion  coefficient  of  species  i  and  /d(=/c  + 
(/jLDe^D^)/(RTctot(De^+x-iD^+x2De^)))  is  the  equivalent  per¬ 
meability  coefficient  of  the  binary  system.  Djjf  and  may  be 
evaluated  by  the  following  equations. 


£  3.198  x  10_8T175 
T  p(v;/3+v73) 
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eff  £  2  I SRT 
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2.3.  Governing  equations  for  electrical  processes 


Electronic  and  ionic  current  densities  are  governed  by  the  charge 
continuity  equations  associated  with  the  Ohm’s  law 


V-iel  =  V-(-afv<pe,) 


V.ii0  =  V.(-crfv^0) 


0 

~Scurrent 
5 current 

S current 

0 


-s 


current 


in  ASL  or  CCCL 
in  AFL 
in  CFL 

in  AFL 

in  electrolyte 
in  CFL 


(7) 


(8) 


where  iet  (iI0)  is  the  electronic  (ionic)  current  density  vector,  (pei 
{(pio)  the  electronic  (ionic)  potential,  SCUrrent  the  current  source.  Lor 
the  charge  conservation,  the  source  of  electronic  current  is  also  the 
sink  of  ionic  current.  oeJ  (<j.°  )  is  the  effective  electronic  (ionic) 
conductivity  of  the  porous  electrode  layer.  The  effective  electronic 
and  ionic  conductivities  of  composite  electrode  are  dependent  on 
structural  parameters  such  as  particle  size,  porosity,  component 
volume  fractions.  The  relationship  between  the  effective  electrode 
properties  and  the  microstructure  parameters  may  be  obtained  by 
the  coordination  number  theory  based  percolation  model  [24-27]. 
Lor  a  binary  system  with  random  packing  of  spheres,  the  effective 
electric  conductivity  of  /c-phase  is  estimated  as  [28], 


<Pk-<Pi  \2 

1  +s/(l  -s) -(/)[) 


0) 


where  is  the  electric  conductivity  of  lt-material  in  dense  solid. 
(j)k  is  the  volume  fraction  of  the  k-particles  in  the  solid  structure, 
the  percolation  threshold  volume  fraction  of  the  k-particles  which 
is  determined  by  [26,29], 


0f//rP/ 
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where  Z  is  the  average  coordination  number  for  each  particle  and 
set  as  6  for  a  random  packing  of  spheres  [29-31  ].  <7°  (£2-1  m-1 )  for 
Ni,  LSM  and  YSZ  are  estimated  respectively  as  [16,32], 


a-”.  =  3.27  x  106-  1065.37 
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The  remaining  quantity  in  Eqs.  (7)  and  (8),  SCUrrent,  may  be  cal¬ 
culated  as, 


where  £  is  the  porosity,  r  the  tortuosity  factor,  rg  the  pore  radius, 
Vi  and  Mj  the  diffusion  volume  and  molar  mass  of  species  i,  respec¬ 
tively.  rei  (rio)  is  the  mean  electronic  (ionic)  conducting  particle’s 
radii,  4>el  (0IO)  the  volume  fraction  of  the  electronic  (ionic)  conduct¬ 
ing  particles  in  the  composite  electrode  material. 

The  molar  mass  conservation  equation  of  species  i  is  given  by: 

f  0  in  ASL  or  CCCL 

V  •  N/  =  <  (6) 

[Ri  in  AFL  or  CFL 

where  Ri  is  the  molar  rate  of  production  or  consumption  of  species 
i  due  to  the  electrochemical  reaction. 


S current  —  ^TPBhrans  (15) 

where  XTPB  is  the  TPB  length  per  unit  volume  (mnrr3)  or  unit  area 
(mm-2)  of  ALL  or  CEL,  and  itrans  is  the  charge  transfer  current  per 
TPB  length.  The  volume-specific  TPB  length  XTPB  can  be  written  as 
[16], 


XjPB  =  27rmin(re/,  ri0)  sin 


^el^el-ioPelPio 


(16) 


where  nvel  =  ((1  -  £)0e/)/(47rr3  /3)  is  the  number  of  electronic  con¬ 
ducting  particles  per  unit  volume  of  the  electrode  function  layer,  6 
is  the  contact  angle  between  particles  (0  is  assigned  as  30°)  [29],  Pk 
( 1<  =  el ,  io)  is  the  percolation  probability  of  k-particle,  Zj_k  (j,  k  =  el, 
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io)  is  the  average  number  of  contacts  between  a  j-particle  and  k- 
particle.  Pk  and  Zj_k  are  calculated  as, 
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The  area-specific  TPB  length,  kjPB,  may  be  written  as 
[16] 


=  4  nrei  sin 


reinvelPei 


(20) 


The  charge  transfer  current  density,  itmns ,  may  be  calculated 
according  to  the  empirical  Butler-Volmer  equation.  For  Ni/YSZ  AFL 
TPBs,  itrans  is  calculated  as  [33,34], 
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For  LSM/YSZ  CFL  TPBs,  itrans  is  calculated  as  [33,34] 
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in  Table  2.  The  contact  resistance  is  set  on  the  interface  between 
the  interconnect  ribs  and  the  electrodes.  That  is,  local  current  den¬ 
sities  cross  the  rib/anode  (iVi^ast)  and  the  cathode/rib  (icca^rib) 
interfaces  are  determined  by 


hib^ASL  = 


( Pe,rib/ASL  ^Pe^SL/rib 

Tcoan 

contact 


(25) 


lCCCL^rib  = 


PeXCCL/rib  ~  Pejib/CCCL 


AZRca 
^  contact 


(26) 


where  ASR“£ntact  ( ASRcc“ntact )  is  the  area  specific  contact  resistance 
at  the  rib-anode  (cathode)  interface,  tpejibiASL  ( VeASLirib )  is  the  elec¬ 
tric  potential  at  the  rib  (anode)  side  of  the  anode-rib  boundary, 
PeXCCLirib  ( Pejib/cccL )  is  the  electric  potential  at  the  cathode  (rib)  side 
of  the  cathode-rib  boundary.  It  should  be  pointed  out  that  the  con¬ 
tact  resistance  is  an  effective  model  parameter  that  may  come  from 
various  origins  such  as  the  oxide  scale  formation  of  the  interconnect 
material  and  the  loose  contact  between  the  cell  components. 


2.5.  Numerical  solution 


The  model  was  implemented  in  the  finite  element  commercial 
software  COMSOL  MULTIPHYSICS®  Version  3.5.  Structured  mesh 
elements  were  used  and  consisted  of  1152  rectangles  with  42,487 
degrees  of  freedom.  The  direct  solver  (UMFPACK)  was  used  to  solve 
the  coupled  partial  differential  equations  of  electronic,  ionic  and  gas 
transport  with  a  relative  convergence  tolerance  of  1  x  10-6.  Subtle 
features  in  the  numerical  solution  process  worthy  mentioning  also 
include:  (1)  As  the  electrochemical  reaction  in  the  region  of  CFL 
underneath  the  rib  may  lead  to  negative  partial  pressure  of  oxygen 
and  cause  numerical  failure,  oxygen  molar  fraction  was  set  as  Xq2  = 
max(xo2 , 10-9)  to  avoid  the  numerical  instability;  (2)  Similarly,  the 
activation  polarizations,  rj™.  and  pcaact,  were  set  to  be  non-negative 
values,  r)™t  =  max(?7^t,  0)  and  rjcaact  =  max(ijcaact,  0),  to  comply  with 
the  physical  requirement. 


Here  a“n  (c^a)  and  {Pcra)  are  the  anode  (cathode)  forward  and 
reverse  reaction  symmetric  factors,  respectively.  Eh2  (Eo2  )  is  the 
activation  energy  for  the  anode  (cathode)  electrochemical  reac¬ 
tion.  Parameters  ij™  and  i™  are  obtained  from  fitting  experiments 
I-V  relations.  rj™t  and  ijcaact  are  respectively  the  anode  and  cathode 
activation  polarizations  and  calculated  as, 

an  an  ET  (  P°h2  Ph2o\ 

nact  =  <Pel  -  (pio  -  Cnc  =  <Pel  ~  <Pio  ~  ln  Xm  (23) 

2F  \P°h2o  Ph2  J 

RT  P^o 

nCact  =  <Pio  -  <Pel  -  ricinc  =  <Pio  ~  (pel  ~  ^  •"(  ^ )  (24) 

where  F  is  the  Faraday  constant,  (p™nc  (^onc)  anode  (cathode) 
concentration  polarization,  p°2  (Ph2o)  e  Partial  pressure  of  H2 

(H20)  at  the  fuel  channel/anode  interface,  the  partial  pressure 
of  02  at  the  air  channel/cathode  interface,  p™  (P™0)  the  Partial 
pressure  of  H2  (H20)  at  the  anode  TPBs,  and  p™  the  partial  pressure 
of  02  at  the  cathode  TPBs. 

2.4.  Boundary  conditions  (BCs) 

As  described  above,  gas  transport  equation,  ionic  conduction 
equation,  electronic  conduction  equation  and  the  Butler-Volmer 
equation  are  taken  into  account  in  this  model.  Proper  settings  of 
the  boundary  conditions  are  required  to  solve  these  coupled  partial 
differential  equations  correctly.  The  boundary  settings  for  the  mass 
transport  equations  are  shown  in  Table  1.  The  boundary  settings 
for  the  electronic  and  ionic  charge  transfer  equations  are  shown 


2.6.  Model  parameters  and  numerical  validation 

The  platinum  or  silver  meshes  are  the  typical  interconnect  for 
single  cell  testing.  The  meshes  may  be  regarded  as  interconnect 
ribs  with  small  pitch  widths.  A  pitch  width  of  0.2  mm  and  a  rib 
width  of  0.05  mm  corresponding  to  the  experimental  description 
were  used  here  for  the  numerical  model  validation.  The  material 
property  parameters  as  deduced  before  [23]  and  the  experimental 
operating  conditions  of  Zhao  and  Virkar  [2]  are  listed  in  Table  3. 
Notice  that,  as  the  fuel  and  air  utilization  rates  are  often  negligible 
in  the  single  cell  testing,  the  same  molar  concentrations  of  gas 
species  are  used  in  the  BCs  of  the  mass  transports  for  obtaining  the 
I-V  relations.  Based  on  the  parameters  in  Table  3  and  assuming  the 
total  contact  resistance  (the  sum  of  ASR™ntact  and  ASR™ntact )  of  the 
cell  are  respectively  0.0132  £2  cm2,  0.0272  £2  cm2  and  0.1  £2  cm2 
for  800  °C,  700  °C  and  600  °C,  the  simulated  I-V  relations  are 
compared  with  the  experimental  results  of  Zhao  and  Virkar  in 
Fig.  2.  As  shown  in  Fig.  2,  the  theoretical  results  agree  with  the 
experimental  measurements  very  well.  It  is  worthy  mentioning 
that  Zhao  and  Virkar  have  estimated  the  total  contact  resistance  to 
be  0.053  £2  cm2  at  800 °C,  which  is  quite  different  from  that  used 
in  our  model.  Nevertheless,  the  parameter  used  in  the  numerical 
model  and  the  experimental  estimate  are  in  fact  consistent  with 
each  other.  The  contact  resistance  estimated  by  Zhao  and  Virkar 
only  corresponds  to  an  equivalent  contact  resistance  for  the  inter¬ 
connect  current  collection  as  they  attributed  the  contact  resistance 
to  be  distributed  over  the  entire  pitch  width  instead  of  only  the  rib 
width.  As  the  current  only  flow  through  the  electrode/rib  interface, 
the  equivalent  contact  resistance  for  our  numerical  model  would 
correspond  to  ASRecZact^(ASRZtact+ASRCc%actKitch/drib  = 
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Table  1 

Boundary  condition  (BC)  settings  for  mass  transports  in  electrodes  (“Insulation”  means  no  flux  through  the  boundary). 


Equations 

Boundary 

ASL/channel  interface 

AFL/electrolyte  interface 

All  others 

BC  type 

H2  molar  concentration 

H20  molar  concentration 

H2  Inward  molar  flux 

H20  Inward  molar  flux 

Fuel  transfer 

BC 

c°  =  p»2/R/T 

C°  =  P°U2o/R/T 

-ian  XA  12 F 

ian  XA  /2F 

trans  IPB' 

Insulation 

Equation 

Boundary 

CCCL/channel  interface 

CFL/electrolyte  interface 

All  others 

BC  type 

02  molar  concentration 

N2  molar  concentration 

02  Inward  molar  flux 

N2  Inward  molar  flux 

Air  transfer 

BC 

c°=p°02/R/T 

c°=P°N2/R/T 

~iLtransXTPB/4F 

0 

Insulation 

Table  2 

Boundary  settings  for  the  electronic  and  ionic  charge  transfer  equations  (“Electric  Insulation”  means  that  the  normal  component  of  the  electric  current  is  zero). 


Equations 

Boundary 
BC  type 

Rib/CCCL  interface 

Reference  potential  (( pe,nb\ccci ) 

Rib/ASL  interface  CFL/electrolyte  interface  AFL/electrolyte  interface 

Reference  potential  (<; Ve,nb\Asi )  Inward  current  flow  Inward  current  flow 

All  others 

Electronic  transfer 

BC 

Vop 

£0  (=1.122V) 

jca  \  A 
tvcms  TPB 

ian  A 
trans  IPB 

Electric  insulation 

Equations 

Boundary 
BC  Type 

Rib/CCCL  interface  Rib/ASL  interface 

CFL/electrolyte  interface 
Interior  current  source 

AFL/electrolyte  interface 
Interior  current  source 

All  others 

Ionic  transfer 

BC 

jca  y  A 

LtransA'TPB 

jan  y  A 
ltransAjTPB 

Electric  insulation 

0.0132  £2  cm2  *  0.2/0.05  %  0.053  £2  cm2,  in  agreement  with 
the  experimental  estimation. 

The  good  agreement  between  the  theoretical  and  experimental 
results  demonstrates  that  the  theoretical  model  is  capable  of  reli¬ 
ably  predicting  the  fuel  cell  performance  in  practical  operations. 
Hence,  the  numerical  model  is  used  in  the  following  to  examine 
the  influence  of  interconnect  rib  on  the  performance  of  planar  SOFC. 
Moreover,  except  for  those  specified  explicitly  for  each  testing  case, 
the  model  parameters  listed  in  Table  3  are  used  throughout  the 
paper. 

3.  Results  and  discussion 

3.2.  Performances  of  cells  with  ribs 

The  influence  of  the  interconnect  configurations  on  the  cell  per¬ 
formance  was  first  examined  with  two  cells  having  different  pitch 


widths  but  the  same  rib  to  pitch  width  ratios.  One  cell  (single  cell) 
has  a  pitch  width  ( dpitch ),  an  anode  rib  width  (cF" )  and  a  cathode 
rib  width  (cFab)  of  0.2,  0.05  and  0.05  mm,  respectively.  The  corre¬ 
sponding  parameters  for  the  other  cell  (stack  cell)  are  2,  0.5  and 
0.5  mm,  respectively.  Even  though  all  the  other  model  parameters 
and  boundary  settings  are  the  same,  as  shown  in  Tables  1-3,  the 
output  current  density  of  the  stack  cell  is  only  6011  Am-2,  a  reduc¬ 
tion  of  19%  from  the  single  cell  output  of  7381  Am-2.  Clearly,  the 
cell  performance  is  strongly  influenced  by  the  pitch  width  and  a 
small  pitch  width  is  beneficial  for  increasing  the  cell  output.  The 
observation  is  consistent  with  the  analytical  and  numerical  analy¬ 
ses  reported  before  [14,17]. 

The  substantial  difference  between  the  performances  of  the  sin¬ 
gle  cell  and  the  stack  cell  may  be  understood  by  the  difference  in 
their  fuel  and  oxidant  species’  distributions.  As  shown  in  Fig.  3a, 
the  distributions  of  H2  and  02  in  the  single  cell  is  reasonably  uni¬ 
form  at  the  electrode-electrolyte  interfaces  as  its  rib  width  (50  p>m) 


Table  3 

Basic  model  parameters  for  the  reference  cell  [23]. 


Operating  parameter 

T 

Vop 

P°h2 

P° 

po2 

p° 

po2 

700  °C 

0.7  V 

0.97  atm 

0.03  atm 

0.21  atm 

Cell  component 

CCCL 

CFL 

Electrolyte 

AFL 

ASL 

Thickness  (|xm) 

50  [2] 

20  [2] 

8  [2] 

20  [2] 

1000  [2] 

e 

45%  [2] 

26%  [2] 

23%  [2] 

48%  [2] 

4>el 

100%  [2] 

47.5%  [2] 

55%  [2] 

55%  [2] 

erf  (SitT1) 

4,05  x  103  (800 °C) 

25.5  (800 °C)  [2] 

1413  (800°C) 

430  (800°C)  [2] 

4.00  x  103  (  700 °C) 

25.0  (700  °C) 

1484 (700°C) 

451  (700°C) 

3.90  x  103  (600 °C) 

24.4  (600  °C) 

1555  (600°C) 

473  (600 °C) 

off  (Sm-1) 

0.36  (800  °C) 

4.2 

0.31  (800  °C) 

0.14  (700  °C) 

(800°C)  [2] 

0.12  (700  °C) 

0.04  (600  °C) 

1.6 

0.04  (600  °C) 

(700°C)  [2] 

0.47  (600 °C) 

rellrio 

1 

1.85 

2.71 

2.71 

r,o  (ixm) 

0.6 

0.2 

0.2 

0.3 

Eo2  (J  mol-1) 

130  x  103  [34] 

af, 

0.65,  0.35 

C,  <A  m"’> 

1.25  x  10-4 

£h2  (J  mol  a) 

120  x  103  [33,34] 

a™,  farn 

1,0.5  [33,34] 

>'S "ref  <A  m~') 

8.0  x  10“3 

T 

3 

3 

3 

3 

ASRcoVt(^m2) 

0.01 

ASRCcontaa  (£2 cm2) 

0.01 
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Fig.  2.  Comparison  of  the  theoretical  and  experimental  I-  V  relationships  at  different 
operating  temperatures.  The  open  symbols  and  solid  lines  denote  respectively  the 
measured  and  the  calculated  results. 


is  comparable  to  the  cathode  layer  thickness  and  very  small  com¬ 
pared  to  the  anode  layer  thickness.  Fig.  3b,  however,  shows  a  quite 
different  situation  for  the  stack  cell  as  its  rib  width  (500  p>m)  is 
over  7  times  the  cathode  layer  thickness.  The  02  concentration  at 
the  cathode  TPB  of  the  stack  cell  decreases  from  2.5  mol  m-3  under 
the  channel  to  0.56  mol  m-3  under  the  middle  of  the  rib.  This  is  an 
important  factor  contributing  to  the  reduced  performance  of  the 
stack  cell  with  a  large  pitch  width.  In  the  other  hand,  the  H2  con¬ 
centration  at  the  anode  TPB  of  the  stack  cell  is  roughly  constant 
at  around  10.7  molnrr3  as  the  rib  width  is  only  about  one  half  of 
the  anode  layer  thickness.  Clearly,  the  influences  of  the  anode  and 
cathode  ribs  on  the  gas  transport  are  characteristically  different.  In 
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Fig.  3.  Distributions  of  hydrogen  and  oxygen  concentrations  (mol  m-3)  in  (a)  single 
cell,  (b)  stack  cell. 


Fig.  4.  The  distributions  of  electric  potentials  in  (a)  single  cell,  (b)  stack  cell. 

an  anode-supported  SOFC,  the  rib  to  pitch  width  ratio  for  the  cath¬ 
ode  should  be  relatively  small  in  order  to  reduce  the  concentration 
polarization. 

Another  factor  contributing  to  the  reduced  performance  of  the 
stack  cell  is  the  increased  ohmic  loss  due  to  the  elongated  current 
conduction  path  associated  with  the  large  pitch  width.  Fig.  4  com¬ 
pares  the  distributions  of  electric  potentials  in  the  single  cell  and 
the  stack  cell.  The  potential  variations  inside  the  cathode  and  the 
anode  of  the  single  cell  are  6  mV  and  19  mV,  respectively.  The  cor¬ 
responding  quantities  for  the  stack  cell  are  respectively  47  mV  and 
40  mV.  The  ohmic  polarization  of  the  stack  cell  anode  is  compa¬ 
rable  to  that  of  the  cathode  and  is  relatively  large.  As  the  effect  of 
rib  width  on  the  anode  concentration  polarization  is  rather  limited, 
there  is  more  room  for  us  to  reduce  the  ohmic  loss  by  increasing 
the  anode  rib  width.  Consequently,  the  desirable  rib  to  pitch  width 
ratio  for  the  anode  should  be  relatively  large  in  comparison  with 
that  for  the  cathode.  Evidently,  the  optimal  rib  to  pitch  width  ratios 
for  the  anode  and  the  cathode  should  not  be  constrained  to  a  single 
value  and  should  be  searched  separately. 

3.2.  The  effect  of  anode  rib  size  on  the  cell  performance 

To  gain  more  insight  into  the  effect  of  the  anode  rib  size  on 
the  cell  performance,  the  stack  cell  output  for  a  fixed  pitch  width 
of  2  mm  is  examined  by  varying  the  rib  width  together  with  a 
few  likely  influential  model  parameters.  The  numerical  results  are 
shown  in  Fig.  5 

The  cell  performance  is  affected  by  both  the  ohmic  and  con¬ 
centration  polarizations.  It  is  natural  to  expect  that  the  contact 
resistance  at  the  rib-ASL  interface  is  an  influential  parameter. 
Many  efforts  are  devoted  to  reduce  ASRcontact  to  the  desired  level 
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Fig.  5.  The  effect  of  anode  rib  size  on  the  cell  output  current  density  (iout)  influenced  by:  (a)  the  contact  resistance,  (b)  the  fuel  composition,  (c)  the  electronic  conductivity 
of  ASL,  (d)  the  porosity  of  ASL,  and  (e)  the  thickness  of  ASL. 


of  0.01  -0.05  Q  cm2  [35].  Consequently,  three  representative  values 
in  this  range,  0.01,  0.025  and  0.05  £2  cm2,  are  used  here.  As  may  be 
seen  in  Fig.  5a,  the  maximal  current  densities  for  ASR“"ntact  of  0.01, 
0.025  and  0.05  £2  cm2  are  6339,  601 8  and  5609  A  nrr2,  respectively, 
corresponding  to  the  optimal  rib  widths  of  1.1,  1.26  and  1.43  mm. 
As  expected,  the  maximal  cell  output  increases  with  the  decrease  of 
the  contact  resistance,  while  the  optimal  rib  width  increases  with 
the  increase  of  the  contact  resistance.  Without  the  rib  size  opti¬ 
mization,  e.g.,  for  a  fixed  rib  width  of  0.4  mm,  the  output  current 
densities  for  ASR“"ntact  of  0.01,  0.025  and  0.05  £2  cm2  are  respec¬ 
tively  7.7%,  14.3%  and  22.7%  less  than  their  optimal  values.  Clearly, 
the  anode  rib  width  has  a  significant  impact  on  the  cell  perfor¬ 
mance.  The  larger  the  contact  resistance  is,  the  more  important  the 
rib  size  optimization  is.  Moreover,  the  benefit  gained  through  rib 
size  optimization  may  substantially  offset  the  shortcoming  associ¬ 
ated  with  the  high  ASRcontact-  For  example,  the  optimized  output  of 
5609  Am-2  for  ASR™ntact  of  0.05  £2  cm2  is  8.7%  higher  than  that  of 
5 160  Am-2  for  ASR“£ntact  of  0.025  £2  cm2  with  a  rib  width  of  0.4  mm. 
This  result  indicates  that  the  advantage  associated  with  the  small 
contact  resistance  may  be  greatly  reduced,  or  even  lost  completely, 
if  the  rib  width  is  not  chosen  appropriately. 


High  fuel  utilization  is  required  in  practical  applications  and 
the  H2  partial  pressure  changes  substantially  along  the  fuel 
channel.  The  rib  size  effect  on  the  cell  performance  is  exam¬ 
ined  here  for  three  representative  H2  partial  pressures,  p°2  = 
0.97,  0.6  and  0.24  atm.  As  shown  in  Fig.  5b,  the  cell  outputs 
with  a  fixed  rib  width  of  0.4  mm  are  7.7%,  9.1%  and  7.0%  less  than 
their  optimal  values  forpf^  of  0.97,  0.6  and  0.24  atm,  respectively. 
The  results  show  that  the  impact  of  the  rib  size  is  not  negligible  for 
a  wide  range  ofpf^.  Nevertheless,  the  influence  ofpf^  on  the  opti¬ 
mal  rib  width  is  moderate.  The  optimal  rib  widths  are  respectively 
1.10,  1 .28  and  1.12  mm  for  p°  of  0.97,  0.6  and  0.24  atm.  The  cur¬ 
rent  output  with  d™b  of  1 .1  mm  for  p°2  =0.6  atm  is  only  0.3%  less 
than  the  optimal  result.  This  is  fortunate  as  the  optimized  anode 
rib  design  is  of  general  applicability  for  different  fuel  operations. 

The  effective  electronic  conductivity  {erf/)  of  ASL  depends 
strongly  on  its  material  composition  and  microstructure  such  as 
the  average  conducting  particle  size  and  volume  fraction  and  the 
porosity.  For  example,  the  measured  electronic  conductivity  of  ASL 
at  800  °C  may  vary  from  430  Snrr1  [2]  to  60774  S  nrr1  [36]  for  seem¬ 
ing  similar  materials.  Therefore,  three  representative  values  of  creJ , 
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Fig.  6.  The  effect  of  cathode  rib  size  on  the  cell  output  current  density  influenced  by:  (a)  the  contact  resistance,  (b)  the  electronic  conductivity  of  CCCL,  (c)  the  porosity  of 
CCCL  and  (d)  the  thickness  of  CCCL. 


45 1 , 45 1 0, 45 1 00  S  rrr 1 ,  are  used  here  to  examine  their  influence  on 
the  rib  size  effect  on  the  cell  performance.  As  shown  in  Fig.  5c,  the 
output  current  densities  for  a  fixed  rib  width  of  0.4  mm  are  5851, 
6503  and  6580Am-2  for  of  of  451,  4510,  45,1005m-1,  respec- 
tively.  With  the  rib  width  optimization,  the  output  may  increase  to 
6339,  6789  and  6846  A  m-2,  or  the  increases  of  8.3%,  4.4%  and  4.0%, 
respectively.  Hence,  the  rib  size  optimization  is  helpful  for  improv¬ 
ing  the  cell  performance  for  any  practical  erff ,  but  the  improvement 
is  more  noticeable  for  low  creJ  of  ASL.  Moreover,  although  the  opti¬ 
mal  rib  widths  for  of  of  451,  4510,  45,100  Snrr1  are  respectively 
1.1,  0.87  and  0.85  mm,  the  rib  width  of  1.1  mm  is  sufficiently  opti¬ 
mal  for  any  oeJ  =451  S  m-1  as  the  output  is  within  1%  of  the 
optimal  value.  Therefore,  the  rib  width  optimized  for  low  oeJ  may 

be  directly  usable  for  the  larger  creef  cases.  Another  point  worthy 
mentioning  is  the  maximal  cell  output  may  increase  by  7%  when 
<reef  is  increased  from  451  to  4510  Snrr1,  but  the  output  increase 

is  less  than  1%  when  aeJ  is  increased  from  4510  to  451 00 Snrr1. 

Consequently,  the  effort  to  increase  aeJ  beyond  4510 Snrr1  is  not 
very  meaningful. 

Fig.  5d  shows  the  output  current  density  as  a  function  of  the 
anode  rib  width  for  ASL  with  three  different  porosities.  The  cell 
outputs  are  seen  to  vary  notably  with  the  rib  width  for  all  the 
three  porosities  examined.  The  optimal  rib  width  increases  with 
the  increase  of  the  porosity  for  reduced  gas  transport  resistance. 
However,  the  change  of  the  optimal  rib  width  with  the  porosity 
is  rather  modest  and  the  optimal  rib  widths  are  1.12,  1.08  and 
1.04  mm  for  the  porosity  of  0.5,  0.45  and  0.4,  respectively.  If  the 
rib  width  optimized  for  one  of  the  three  porosities  is  used  for 
an  ASL  with  the  other  porosity,  the  current  output  is  different 
from  its  optimal  value  by  less  than  0.5%.  Therefore,  the  optimal 
rib  width  may  be  viewed  as  independent  of  the  ASL  porosity  in 
practice. 


Fig.  5e  shows  the  relationship  between  the  output  current  den¬ 
sity  and  the  anode  rib  width  for  three  ASL  thicknesses.  The  rib  size 
effect  on  the  cell  output  is  significant  for  all  the  three  ASL  thick¬ 
nesses.  The  ASL  thickness  also  has  a  significant  effect  on  the  cell 
output.  However,  the  optimal  rib  widths  for  the  three  thicknesses 
are  roughly  the  same,  agreeing  with  the  previous  observation  that 
the  optimal  rib  width  is  independent  of  the  ASL  thickness  as  long 
as  the  thickness  is  larger  than  a  few  hundred  microns  [19]. 

3.3.  The  effect  of  cathode  rib  size  on  the  cell  performance 

Fig.  6a  shows  the  results  on  the  variation  of  the  output  current 
density  with  the  cathode  rib  width  for  different  contact  resistances. 
The  maximal  current  densities  for  ASRcc“ntact  of  0.01,  0.025  and 
0.05  £2  cm2  are  6035,  5425  and  4785  Am-2,  respectively,  corre¬ 
sponding  to  the  optimal  rib  widths  of  0.46,  0.56  and  0.67  mm.  As 
discussed  above,  the  optimal  rib  width  for  the  anode  is  in  the  vicin¬ 
ity  of  1  mm.  Compared  to  the  results  obtained  with  a  fixed  cathode 
rib  width  of  1  mm,  the  rib  width  optimization  increases  the  out¬ 
put  current  by  20%,  13%  and  7.1%  for  ASRcc“ntact  of  0.01,  0.025  and 
0.05  £2  cm2,  respectively.  Clearly,  the  cathode  rib  width  has  a  sig¬ 
nificant  impact  on  the  cell  performance.  Moreover,  the  optimal  rib 
widths  for  the  cathode  and  the  anode  are  quite  different. 

Fig.  6b  shows  the  results  of  the  output  current  densities  for 
three  representative  electronic  conductivities  of  CCCL,  creJ  =  800, 
4000,  20,000 Snrr1.  The  output  current  densities  for  a  fixed  cath¬ 
ode  rib  width  of  1  mm  are  4370,  5036  and  5289  Am-2  for  of 
800, 4000, 20,000  S  nrr 1 ,  respectively.  With  the  rib  width  optimiza¬ 
tion,  the  outputs  may  increase  to  4895, 6035  and  6453  A  nr  2 ,  or  the 
increases  of  12%,  20%  and  22%,  respectively.  Evidently,  the  cathode 
rib  size  optimization  is  very  helpful  for  improving  the  cell  perfor¬ 
mance  for  any  practical  aeJ ,  and  the  benefit  often  increases  with  the 

increased  aeJ  of  CCCL.  As  shown  in  Fig.  6b,  however,  although  the 
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optimal  cell  output  increases  with  the  increased  aeJ  of  CCCL,  the 

optimal  cathode  rib  widths  for  aeJ  of  800, 4000,  20,000  S  nrr 1  only 
vary  slightly  and  are  respectively  0.49,  0.46  and  0.45  mm.  In  fact, 
the  rib  width  of  0.46  mm  is  sufficiently  optimal  for  all  aeJ  exam¬ 
ined  as  the  output  is  within  0.1%  of  the  optimal  value.  Therefore, 
the  rib  width  optimized  for  a  medium  creJf  of  CCCL  (~4000Sm_1) 

is  representative  for  all  practical  aeJ  cases. 

Fig.  6c  shows  the  output  current  density  as  a  function  of  the 
cathode  rib  width  for  CCCL  with  three  different  porosities.  Similar 
to  the  case  for  the  optimal  anode  rib  sizes,  the  change  of  the  optimal 
cathode  rib  width  with  the  porosity  is  rather  modest  and  the  opti¬ 
mal  cathode  rib  widths  are  0.44,  0.46  and  0.48  mm  for  the  porosity 
of  0.4,  0.45  and  0.5,  respectively.  If  the  rib  width  optimized  for  a 
porosity  of  0.4  is  used  for  a  CCCL  with  a  porosity  of  0.5,  the  current 
output  is  different  from  its  optimal  value  by  less  than  1%.  There¬ 
fore,  the  optimal  cathode  rib  width  may  be  viewed  as  independent 
of  the  CCCL  porosity. 

Fig.  6d  shows  the  relationship  between  the  output  current  den¬ 
sities  and  the  cathode  rib  width  for  three  different  CCCL  thicknesses. 
As  shown  in  the  figure,  the  rib  size  effect  on  the  cell  output  is  signif¬ 
icant  for  all  the  three  CCCL  thicknesses.  The  CCCL  thickness  also  has 
a  significant  effect  on  the  cell  output.  The  maximal  cell  outputs  are 
6035,  6399  and  6528  Am-2  for  the  CCCL  thicknesses  of  0.05,  0.10 
and  0.15  mm  and  are  20%,  15%  and  11%  higher  than  the  outputs 
with  a  fixed  cathode  rib  width  of  1  mm,  respectively.  Some  depen¬ 
dence  of  the  optimal  cathode  rib  width  on  the  CCCL  thickness  is 
observed.  The  optimal  rib  widths  are  0.46,  0.55  and  0.6  mm  for  the 
CCCL  thickness  of  50,  100  and  150  |jim,  respectively.  However,  the 
dependence  of  the  optimal  cathode  rib  width  on  the  CCCL  thick¬ 
ness  is  rather  mild.  If  the  rib  width  optimized  for  a  CCCL  thickness 
of  50  p>m  is  used  for  a  CCCL  thickness  of  150  p>m,  the  current  out¬ 
put  is  different  from  its  optimal  value  by  less  than  1%.  Therefore, 
the  cathode  rib  width  optimized  for  the  CCCL  thickness  of  50  p^m  is 
sufficiently  optimal  for  other  practical  CCCL  thickness. 
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Fig.  7.  Dependence  of  the  optimal  anode  (cathode)  rib  width,  driban  (dr;faca),  on  the 
contact  resistance  and  pitch  width:  (a)  anode,  (b)  cathode. 


As  described  above,  Eq.  (30)  is  only  an  approximation  for  the  true 
numerical  results.  The  usefulness  of  Eq.  (30)  for  guiding  the  optimal 
anode  rib  width  design  may  be  vilified  by  the  accuracy  parameter, 
y,  defined  as, 


3.4.  Expressions  for  the  optimal  rib  widths 


Y  = 


Imax  ~  bib 
l  max 


x  100% 


(31) 


As  discussed  above,  the  output  current  density  depends  strongly 
on  the  rib  width  and  a  suitable  choice  of  the  rib  width  is  very  impor¬ 
tant  for  the  high  performance  of  a  SOFC  stack  cell.  The  numerical 
results  shown  above  also  demonstrate  that  the  optimal  rib  width 
for  a  given  pitch  width  is  only  sensitive  to  the  rib-electrode  con¬ 
tact  resistance.  That  allows  us  to  obtain  simple  expressions  for  the 
optimal  rib  widths,  providing  an  easy  to  use  guidance  for  the  broad 
SOFC  engineering  community. 

The  optimal  rib  width  for  a  given  contact  resistance  and  a  given 
pitch  width  may  be  obtained  by  changing  the  rib  width  to  achieve 
the  maximum  cell  current  density.  The  optimal  rib  widths  for  prac¬ 
tical  contact  resistances  and  pitch  widths  are  shown  in  Fig.  7.  As 
may  be  seen  from  Fig.  7a,  the  optimal  anode  rib  width  for  a  given 
pitch  width  is  dependent  on  the  contact  resistance  approximately 
linearly.  Moreover,  the  slope  in  the  linear  relationship  is  roughly 
independent  of  the  pitch  width.  Similarly,  the  optimal  anode  rib 
width  for  a  given  contact  resistance  is  also  found  to  be  dependent 
approximately  linearly  on  the  pitch  width  with  the  slope  indepen¬ 
dent  of  the  contact  resistance.  Consequently,  the  optimal  anode  rib 
width  may  be  expressed  as: 

drib  =  A  +  B  x  dpjfth  +  C  x  ASRContact  (30) 

The  three  parameters  in  Eq.  (30)  for  the  optimal  anode  rib 
width  may  be  obtained  by  fitting  as  Aan  =  0.21mm,  Ban  =  0.42, 
Can  =  10  mm  cm-2.  As  will  be  discussed  below,  Eq.  (30)  is  also 
suitable  for  estimating  the  optimal  cathode  rib  width,  but  the 
three  parameters  should  be  revised  as  Aca  =  0.13  mm,  Bca  =  0.14,  and 
Cca  =  5.5  mm  cm-2. 


where  imax  is  the  output  current  density  corresponding  to  the 
numerically  optimized  rib  width,  while  irib  is  the  current  density 
for  the  rib  width  determined  analytically  by  Eq.  (30).  Fig.  8a  shows 
the  results  of  y  for  different  anode  pitch  widths  and  contact  resis¬ 
tances.  Clearly,  y  is  always  negligibly  small.  Therefore,  the  anode 
rib  width  determined  by  Eq.  (30)  is  highly  optimal  for  engineering 
applications. 

As  shown  in  Fig.  7b,  the  optimal  cathode  rib  width  for  a  given 
pitch  width  is  dependent  on  the  contact  resistance  approximately 
linearly,  but  the  slope  increases  with  the  pitch  width.  Similarly, 
the  optimal  cathode  rib  width  for  a  given  contact  resistance  is 
also  found  to  be  dependent  approximately  linearly  on  the  pitch 
width  and  the  slope  increases  with  the  contact  resistance.  It  may 
be  inferred  then  that  Eq.  (30)  may  not  accurately  represent  the  opti¬ 
mal  cathode  rib  width.  With  the  rib  and  pitch  widths  in  mm  and 
ASR  in  ficm2,  the  cathode  rib  width  is  best  fitted  as, 

C  =  0.28  +  0.07  x  dc“tch  +  (0.8  +  2.35  x  dcpaitch)  x  ASRccaontact  (32) 

However,  considering  the  extremal  property  of  the  cell  output 
around  the  optimal  rib  width,  the  deviation  of  the  cell  output  from 
its  optimal  value  is  only  quadratic  on  the  small  error  of  estimating 
the  optimal  rib  width.  It  is  likely  then  that  the  estimation  of  the 
optimal  rib  width  by  the  simple  expression  of  Eq.  (30)  with  proper 
parameters  suffices  for  any  practical  purpose.  To  this  end,  the 
fittings  yield  Aca  =  0.13  mm,  Bca  =  0.14  and  Cca  =  5.5  mm  cm-2. 
Fig.  8b  shows  the  deviation  of  the  current  output  with  the  optimal 
cathode  rib  width  estimated  by  Eq.  (30)  from  that  with  truly  opti¬ 
mized  rib  width  for  various  pitch  widths  and  contact  resistances. 
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Fig.  8.  Deviations  of  the  current  outputs  with  analytically  determined  rib  widths 
from  that  with  numerically  optimized  rib  widths  for  different  pitch  widths  and 
contact  resistances:  (a)  anode,  (b)  cathode. 


As  shown  in  Fig.  8b,  the  maximal  error  by  using  Eq.  (30)  is  less  than 
1%  for  all  cases  examined.  The  result  demonstrates  that  the  optimal 
cathode  rib  width  may  also  be  accurately  estimated  by  Eq.  (30). 

4.  Summary 

We  have  described  a  comprehensive  two  dimensional  math¬ 
ematical  model  for  the  performance  of  the  planar  SOFC  stack. 
The  model  takes  into  account  the  contact  resistances  between  the 
electrodes  and  the  interconnect  ribs  and  the  dependence  of  the 
effective  electrode  properties  on  the  microstructure  parameters 
of  the  porous  electrodes.  The  impacts  of  the  electrode  rib  widths 
are  systematically  examined  by  varying  the  contact  resistance,  fuel 
composition,  electrode  porosity,  thickness  and  conductivity.  The 
numerical  results  show  conclusively  that  the  output  current  den¬ 
sity  depends  strongly  on  the  rib  widths  and  the  optimal  rib  widths 
for  the  anode  and  the  cathode  are  quite  different.  Nevertheless, 
the  optimal  widths  for  both  the  anode  and  cathode  ribs  are  found 
to  be  only  sensitive  to  the  contact  resistance  and  the  pitch  width 
and  may  be  determined  by  the  same  simple  analytical  expres¬ 
sion,  =A  +  B  x  dp itch  +  C  x  ASRContact •  For  the  anode,  A  =  0.21  mm, 
B  =  0.42,  C=  1 0  mm  cm-2  and  drib  is  the  optimal  anode  rib  width 
corresponding  to  a  given  anode  pitch  width,  dpitch,  and  a  contact 
resistance  at  the  interconnect  rib  and  anode  interface,  ASRC0ntCict- 


For  the  cathode,  A  =  0.13  mm,  B  =  0.14,  C=  5.5  mm  cm-2  and  drib 
is  the  optimal  cathode  rib  width  corresponding  to  a  given  cathode 
pitch  width,  dpitch,  and  a  contact  resistance  at  the  interconnect  rib 
and  cathode  interface,  ASRcontact •  The  simple  formula  provides  an 
easy  to  use  guidance  for  the  broad  SOFC  engineering  community  in 
designing  the  optimal  rib-channel  layout. 
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